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i— K Abstract 

■ A model of inhomogeneous baryogenesis based on the Affleck and 

■ Dine mechanism is described. A simple coupling of the scalar baryon 
field to the inflaton allows for formation of astronomically significant 

£^ ■ bubbles with a large baryon (or antibaryon) asymmetry. During the 

farther evolution these domains form compact stellar-like objects, or 
lower density clouds, or primordial black holes of different size. Ac- 

OO ! cording to the scenario, such high baryonic number objects occupy 

relatively small fraction of space but despite that they may signifi- 
£> . cantly contribute to the cosmological mass density. For some values 

of parameters the model allows the possibility the whole dark matter 
S_i ■ in the universe to be baryonic. Furthermore, the model allows the 

existence of the antibaryonic B-bubbles, i.e. a significant fraction of 
the mass density in the universe can be in the form of the compact 
antimatter objects (e.g. anti-stars). 



1 Introduction 

Though the universe in our neighbourhood looks as 100% baryonic, the ob- 
servations do not exclude sub dominant but astronomically large domains and 
objects of antimatter. For a recent analysis and literature see e.g. ref. PQI2]- 
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The standard mechanism of baryogenesis [3] predicts baryon asymmetry ex- 
pressed as the ratio of baryon-to-photon number densities: 



P = ub/ti. 
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(1) 



to be constant both in space and, after baryogenesis is over, also in time, 
however, more complicated models can lead to not so dull pattern. The 
outcome of the baryogenesis scenario is closely related to the mechanism of 
CP- violation realized in cosmology For a review of the latter see ref. [I]. 
If CP is broken spontaneously [5], the universe would be charge symmetric 
with equal weight of baryonic and antibaryonic domains so that the average 
value of (3 over whole universe would be zero [61 [7] . However, in the standard 
cosmology the domain size would be very small by astronomical standards 
and some exponential (inflationary) expansion is necessary to avoid immedi- 
ate problems [HI E] . Another cosmological problem of domain walls with huge 
energy density inside our horizon [6] can be solved either by non-restoration 
of CP invariance at high temperature [10J (TTJ, but in this case antimatter 
domains would not be created, or by a possible mechanism of wall destruc- 
tion [T2l 13J. However, even if these problems are resolved, it can be shown 
that baryo-symmetric cosmology demands the domain size above (a few) 
Gpc. Otherwise too intensive gamma ray background would be created by 
the annihilation on the domain boundaries [H] . 

If in addition to spontaneous CP- violation there exists an explicit one, the 
universe would be baryo-asymmetric and, if the magnitude of spontaneous 
CP breaking happened to be larger than the explicit one, there could exist 
subdominant antimatter domains p21 EH EE] • 

Here we explore another mechanism of inhomogeneous baryogenesis sug- 
gested in ref. [TB]. This mechanism is based on the Affleck and Dine [T7j 
idea that in supersymmetric theories a complex scalar field x with non-zero 
baryonic number may be displaced far from the origin, \ = 0, along flat 
directions of the potential U (%) . Such flat directions are generic features of 
supersymmetric models. Inflationary stage is especially favorable for such a 
travel along the flat directions due to the infrared instability of massless or 
light, m 2 , <C H 2 , fields in De Sitter space-time. According to the results of 
ref. [THl EU EH] ( see however, ref. [21]), vacuum expectation value of massless 
scalar field rises as: 



Inflation not only drives x away from the origin but also allows for formation 
of the exponentially big domains with a large value of x- 

After the Hubble parameter, H, dropped below m x , the scalar baryon 
(AD-field) x would evolve down to x = transforming the accumulated 



(X 2 ) = Hh/^ . 



(2) 



2 



baryonic charge into quarks by 5-conserving decays. Typically such scenario 
leads to a large baryon asymmetry, j3 ~ 1 and theoretical efforts should be 
taken to diminish it down to the observed value (3 ~ 5 ■ 1CT 10 . 

The modification of this mechanism, according to ref. |16| . was an intro- 
duction of the quadratic coupling of x to the inflaton field $ in the following 
renormalizable form, Ai|x| 2 ($ 2 + A^&), where ji is a constant parameter with 
dimension of mass. This term can be rewritten as: 

^( X ,$) = A 1 |x| 2 ($-$i) 2 (3) 

by redefinition of the mass term of x- 

When the inflaton field evolves down to zero, it passes through $ = $i 
and at that moment the window to flat direction in x potential would be open 
but only for a rather short period of time. Correspondingly, the probability 
for x to evolve away from the origin would be small and baryogenesis resulting 
in high (3 could proceed only in a minor fraction of space. So we should 
expect that in the bulk of space the usual homogeneous baryogenesis took 
place, and compact B-bubbles, with very high baryonic number density could 
be created over this background. Despite their small size, B-bubbles might 
contribute significant amount into the total cosmological mass density and 
even make all dark matter of the universe. They would either form compact 
stellar- like objects, or high density clouds, or primordial black holes (BH). 
In the simplest version of the model there should be an equal number of 
baryonic and antibaryonic B-bubbles. So cosmology would be almost baryo- 
symmetric but the bounds derived in ref. [H] are not applicable here because 
of compactness of the bubbles. An interesting feature of this scenario is 
that dark matter (DM) can be made from the usual particles, baryons and 
antibaryons in the form of evolved stellar type objects or primordial black 
holes. 

The paper is organized as follows. In the Sec. [2] we give a detailed 
description of the model and it's parameter space. The B-bubble evolution 
and the generation of the baryon asymmetry are discussed in Sec. [3J The next 
section, Sec. HJ gives a detailed description of the development of quantum 
fluctuations of the Affleck and Dine field and the distribution of B-bubbles 
over the length and mass. In Sec. [5] the fate of B-bubbles during different 
cosmological epochs is considered. In SecOthe observational bounds on the 
dark matter in the form of primordial black holes (PBH) and the amount af 
antimatter in the compact stellar objects are presented. Our conclusions are 
hosted in the Sec. UJ 
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2 The model and the parameter space 



First we will study qualitatively the evolution of x an d formation of the 
bubbles with a large value of this field. Though the model and results are 
similar to those of ref . [16] , here we discuss a somewhat simpler scenario. 

Let us consider a toy model which can lead to the features described 
in the Introduction. We assume that the evolution of the inflaton field is 
governed by the usual potential: 

= m|$ 2 /2 + A*$ 4 / 4 > (4) 
while x "lives" in the potential 

U X ( X , *) = A x ($ - $!) 2 |x| 2 + A 2 | X | 4 ln K + ml\ X \ 2 + m\ X 2 + m? X * 2 . (5) 

The first term in this expression describes the interaction between the inflaton 
and x fields introduced in ref. [TS]. We assume that Ai and |%| are sufficiently 
small, and so this term has negligible impact on the inflaton evolution. The 
second term is the Coleman- Weinberg potential [22] obtained by summation 
of one loop corrections to the quartic potential, A 2 |x| 4 - The l &s t two mass 
terms are not invariant with respect to the phase rotation: 

X - e l9 X (6) 

and thus break baryonic charge conservation (see below). If Wq < 2|m!| 2 , 
the potential created by these two terms is not bounded from below but the 
effective mass: 

m 2 eff = Ai($ - $x) 2 + m 2 + 2|mi| 2 cos(2a + 26) (7) 

is almost always positive, except for a short time when $ w <&x. Anyhow 
the A2|x| 4 hr|x| 2 t erm stabilizes the potential but its minimum shifts from 
X^ = to some non-zero x- m the above equation a and 9 are the phases 
of mass, m 1 = \mi \ exp(ia), and of AD- field, x — \x\ exp(z6 l ), respectively. 

It is interesting that potential (jSJ) conserves CP despite presence of the 
complex mass parameters. Indeed the imaginary part of m can be excluded 
by the redefinition x ~^ X ex P( — ^ a )- However, a coupling of x to fermions 
may destroy CP-conservation. 

Classical evolution of homogeneous fields $ and x is determined by the 
equations 

$ + 3H$ + m|$ + A$$ 3 + 2Ai|x| 2 (^ - $i) = , (8) 
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X + SHx + Ax($ - $!) 2 X + A 2 | X | 2 X (2 In ^ + lj + 2m*Y = , (9) 

which we solved numerically, but qualitative behavior of the solutions can 
be understood from the following simple considerations. We start from in- 
flationary stage when the cosmological energy density is dominated by the 
inflaton field. We assume that the product Ai|x| 2 is sufficiently small, so the 
X-dependent term in eq. (JHJ) is not essential. In this regime the evolution 
of $ is simply a slow roll in the corresponding potential (j3J) till the Hubble 
parameter remains large in comparison with the potential slope: 

In the course of the cosmological expansion $ approaches the equilibrium 
point $ = and H drops down so that 

H 2 < m 2 ($) =m% + 3A$$ 2 . (11) 

Inflation stops and $ starts to oscillate around $ = producing particles and 
heating the universe. Correspondingly, the De Sitter expansion regime with 
domination of the vacuum-like energy turns into the matter domination and 
later into the radiation domination. To be more precise the equilibrium value 
of $ depends upon x- But as we will see below, x ultimately tends to zero 
and during inflationary period this correction is not of much importance. If 
the life-time of x is longer than the life-time of the inflaton, the impact of x 
on the inflaton behavior might be essential at the (re/pre)heating. 

According to our assumption and with our choice of the parameters, the 
impact of x on the inflaton behavior is almost unimportant. Therefore during 
the qualitative discussion we consider the evolution of x i n the fixed (but 
still time dependent) inflaton background. However, we have not neglected 
corresponding term when solved the equations numerically. 

At the beginning of inflation, $ is large, $ > mpi, and the effective mass 
of X © is also large. During this period potential (jSJ) has only one minimum 
at x — and if m e ff, eq. ((7j), is large in comparison with H, quantum 
fluctuations of x are n °t essential and field x "sits" in the minimum x^ — 0. 

In fact the picture is more complicated. If ^ 0, then for large <3> the 
Hubble parameter would be larger than m e ff and x could quantum fluctuate 
away from zero climbing up the potential slope during an early inflationary 
stage. With decreasing H field x would return back to zero if the classical 
rolling down of x is faster than the decrease of H . 

When the amplitude of the inflaton drops down and approaches $i, the 
effective mass of x decreases and the potential gets the second minimum at 
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some 7^ 0- This second minimum appears and disappears when vr? e tc 
passes down and up the value 2A2<T 2 /e 3 / 2 . Initially the second minimum at 
X = X^ lays above the minimum at x = 0> and later with the diminishing 
difference |<& — $i| it becomes the absolute minimum, but separated from 
the minimum at x = by a barrier. In this situation the first order phase 
transition by quantum tunneling from x — to x — X^ is possible, but the 
transition probability is exponentially small and the process can be neglected. 
The transition to % = x would be much more efficient when the barrier 
became low or disappeared. It is essential that the low barrier existed only for 
some finite time. After $ passed $i and traveled further to zero, the potential 
barrier would reappear. The behavior of U(x, is depicted in Fig. [TJ As 
a result of such non-monotonic behavior of the potential, x would be either 
forced back to the first minimum at x — or would remain beyond the barrier 
near x = X ■ The separation depends upon the amplitude which x managed 
to reach to the moment of the "gate closing". The boundary/critical value 
of x is determined by the position of the maximum of U(x, 3>) at x — X^ 3 \ 
which separates the two minima (see Fig. [1]). x^ is the smaller solution of 
the equation: 

2A 2 | X rin^! + m e 2 // = (12) 

and by an order of magnitude: x^ ~ m eff/\^2- An important feature of 
X3 is that it is not constant but moves together with $. 

Thus we expect that for some time, when the second minimum at x^ 
existed, there should be two phases with large x ~ X and small x ~ 0. 
Which phase dominates is determined by the duration of the "open gate" 
period and the efficiency of quantum fluctuations of x- We will show in 
what follows (Sec. H]) that one can naturally realize the scenario when the 
probability of formation of large x bubbles is sufficiently small and they 
would occupy a small fraction of the total volume of the universe. Still their 
size can be astronomically large if they were formed during inflationary stage. 
In other words, when $ was close to $i, inflation lasted sufficiently long to 
stretch up the bubble size exponentially. In the future such bubbles would 
form objects (or domains) with a large value of the baryon asymmetry and 
we will call them B-bubbles or just BBs. 

Note that the phase transition (p.t.) which we discuss here is neither 
first order, nor second order but something in between. It takes place when 
the barrier is absent and in this sense it is second order but it resulted 
in formation of bubbles of new phase inside the old one and, as such, it 
resembles the first order p.t. Note that the evolution of x demonstrates the 
hysteresis behavior. Initially when x was near zero, it remained there despite 
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Figure 1: Qualitative behaviour of U x {\) for different values of m 2 ,Jt). 



appearance of the second minimum of the potential at x • At a later stage 
when the second minimum became higher than the first one, \ remained 
there and rolled down to x = only when the minimum x^ completely 
disappeared. 

We assume that wig > 2m 2 , so the local minimum of U x never disappears 
and classically x should always remain zero. However, quantum fluctuations 
of light, m < H, scalars in De Sitter background allow to overcome this 
barrier. Roughly speaking it is possible when Ai($ — $i) 2 < H 2 and thus 
the characteristic time when the gate to the negative slope of the potential 
was open, lasted approximately 

^ 3# 2 = 27^(2774 + A»tf) 

~ v/A^i y/%m 2 Pl (m% + A„$ 2 ) ' 1 ' 

where we used the slow roll approximation: <E> ps U$/3H. According to eq. 
(j2J), during this time interval x can reach the value of the order of 



, H 3 At PhT $f (2m| + A$$ 2 ) 5 / 2 



X 4vr 2 V 27Ai m| + A$$ 2 1 J 

For the chosen below parameters y_ can reach the boundary value Xc after 
which it would evolve classically down to the second minimum of U x . 



Since during this period m 2 was much smaller than H 2 , the phase of 
X equally populated the interval [0, 2ir]. This would allow for generation 
of bubbles with a large baryon asymmetry. On the other hand, if x stays 
in the minimum over the phase variable i.e. at cos(2a + 26) = — 1, the 
baryon asymmetry would be small, because no "angular momentum" can 
be generated in the process of the evolution of x down to zero, despite the 
non-sphericity of the potential (see Sec. [3]). Such could be realized if 

X evolved to X2 along the path a + 9 = ir/2. However, even in this case the 
X-state initially localized near this value of the phase, would be dispersed 
by quantum fluctuations over all values of the phase because the mass of 
the pseudogoldstone degree of freedom is much smaller than the Hubble 
parameter. 

In flat space-time the tunneling from x^ = to a nonzero x i n the 
potential A|x| 4 is determined by the Fubini instanton [23], which leads to a 
very low probability of the transition, ~ exp(— 327r 2 /3A). The exponential 
suppression for small A is related to a non-zero energy density of the bubble 
wall, where x interpolated from x — to x — X ■ This suppression is absent 
in our case by the following reasons. First, the effective coupling strength, 
A2 ln(|x| 2 /cr 2 ), in our case is not constant but diverges at x = 0. Second, 
the exponential expansion stimulates rise of quantum fluctuations and leads 
to a large deviations from x — 0. Moreover, all space inhomogeneous effects 
are exponentially smoothed down by the cosmological expansion. A detailed 
discussion of tunneling in curved space-time can be found in ref. [21] 

Another danger comes from the opposite direction. The probability of BB 
formation may be too large and practically all the space would be occupied by 
the large x bubbles, while low x values would be practically absent. However, 
this problem can be easily cured by an adjustment of the parameters of the 
model, i.e. of the masses and coupling constants in potentials (jl]) and 
A sufficiently large term m 2 ,^ 2 can lead to the necessary suppression of the 
bubble formation. 

The picture that we described above can be realized with the following 
rather natural values of the parameters. Since we want inflation to continue 
when $ passes $1, we need $1 > mpf, $ x = (2— 3)mpi is sufficient for creation 
of astronomically large B-bubbles. The inflaton mass should be bounded 
by 10 -7 mpi < m$ < lCT 6 mp; to avoid too large density perturbations, 
Sp/p < 10~ 5 [25]. The coupling constant A$ must be in the interval 10~ 14 < 
A$ < 1CT 12 by the same reason to suppress the density perturbations which 
at horizon crossing are, roughly speaking, equal to Sp/p ~ 10 2 ^/\$ [25] . The 
coupling of x f° the inflaton, Ai, is bounded by Ai < 1CT 6 , because otherwise 
the one loop correction to the inflaton potential induced by this coupling 
would give rise to unacceptably large A$ ~ A 2 /47r 2 . Interesting results are 
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obtained with Ai ~ 10 -10 . Such a small value of the inflaton coupling to 
X naturally fits the idea of very weak coupling of the inflaton to the usual 
matter. We choose the x-self-coupling constant A2 to be near 1CT 2 to obtain 
cosmologically interesting consequences. 

3 B-bubble evolution and baryogenesis in the 
very early universe 

When $ dropped noticeably below $1, the potential barrier reappeared, the 
formation of the bubbles with x ~ X^ from x = terminated, and the 
primeval cosmological matter would consist of the inflaton field, which still 
continued to inflate the universe, with exponentially expanding bubbles of 
non-zero x over practically homogeneous background with x — 0. The energy 
density of BBs were initially negligible with respect to the total cosmological 
energy density, dominated by $, and in a good approximation the universe 
can be considered as homogeneous. 

Further decrease of the inflaton field leads to a large rise of the effective 
mass of x, eq. ([7]), so the second minimum of the potential U (%, $) disappears 
(see Fig. [T]) and % started to roll down to zero. The evolution of the spatially 
homogeneous field x = xif) an d generation of baryon asymmetry can be 
easily visualized with the following mechanical analogy. The equation of 
motion for \ m arbitrary potential reads: 

X + 3H X + U'(x) = 0- (15) 

It is exactly the equation of motion of a point-like particle in potential U(x) 
in Newtonian mechanics. The second term originating from the cosmologi- 
cal expansion corresponds to the liquid friction force. In this language the 
baryonic charge density is simply the angular momentum of the particle 

= i(x*d^ X - 9,X*X) = -2\x\ 2 d,9 . (16) 

If the potential is rotationally symmetric in the complex x plane, the angular 
momentum is evidently conserved. This is the reason for the introduction 
of the rotation non-invariant mass terms in the potential: they break the 
baryonic current conservation: 

= 2i(m?x* 2 - mix 2 ) = 2\ X \ 2 \m x \ 2 sin2(# + a) ^ 0, (17) 

where 9 and a are the phases of x an d wii, respectively. 

Too large values of mi may change the essential features of the potential 
described in the previous section, as it is shown on the Fig. [2j In particular, 
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the negative term in m\jt ((7j) should be smaller than Ai$f to forbid traveling 
of x from zero to larger values when $ is far from $ x . A stronger bound on mi 
follows from the requirement that the classical slow roll along the minimum 
at 2 (a + 6) = it, is not too fast. According to eq. ([9]) this would be avoided 
if mi is bounded from above by \m\\ < m$. The classical roll down from 
X = would be absent for any $ if m^ > 2|mi| 2 




Figure 2: The form of U x in the complex %-plane near the origin. On the 
left: mi is small and the potential is almast symmetric. On the right: mi 
is big and this changes the form of the potential, the minimum becomes a 
saddle point. 

For such values of the parameters the non-sphericity of the potential is 
quite small and there is no preferred direction in the complex %-plane for 
the quantum fluctuations generated at inflationary stage when m e ff < H. 
On the other hand, when inflation is over and x approaches its equilibrium 
point at x = 0, the non-sphericity of the potential is essential for creation 
of the noticeable angular momentum, i.e. for generation of baryonic charge 
asymmetry in the sector of the scalar baryon x, see Fig. [31 

Thus to summarize. As it is already described above, during inflation, 
when $ is close to $i, the quantum fluctuations of x become large enough 
to overshoot the barrier and field x starts classical rolling down towards the 
deeper minimum, x {t)- As one can see m Figs. [3] and HJ x reaches this 
minimum, follows its evolution in time, and oscillates around this minimum 
with decreasing amplitude. At the end of inflation when the minimum of the 
potential at x = becomes again the only minimum, the field rolls down to 
the origin along the slope of the potential. Coming close to the origin the 
field starts to rotate around it due to the angular momentum generated by 
non-sphericity of the potential, i.e. due to the complex mass terms (see Fig. 
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Figure 3: The typical picture of the evolution of x m the complex %-plane, 
X = Xi + iX2- The zoomed part inside the figure shows the rotation around 
the absolute minimum after the end of inflation. 

[3]and the zoomed part inside). To gain a considerable angular momentum 
the magnitude of this terms at low x should be comparable with the others 
entering equation of motion (J9j). This is the reason why we have chosen mo 
only slightly larger than \mi\. The term ~ x 4 is not essential near x — 0. The 
form and slope of the spiral-like trajectory depend on the effective mass of x- 
The rotation of x around the origin corresponding to the nonzero baryonic 
charge would be transformed into the baryonic asymmetry of quarks and 
ultimately of protons and neutrons by B-conserving decays of x into (light 
or massless) fermions. The sign of the asymmetry depends on the direction 
of the rotation and can be arbitrary because, as we mentioned above the 
initial value of the phase 9 is uniformly distributed in the interval [0, 2tt]. So 
in this version of the scenario we would expect an equal number of baryonic 
and antibaryonic bubbles. In slightly complicated models this symmetry 
between B and anti-B bubbles may be lifted. 

The magnitude of the asymmetry inside the B-bubbles can be much larger 
than canonical value ([T]). Moreover, if the life-time of x noticeably exceeds 
the life-time of the inflaton, the asymmetry may be even larger than unity, 
because the energy density of the decay products of the inflaton into rela- 
tivistic particles would be strongly red-shifted with respect to non-relativistic 
longer-lived x-bosons. The baryon asymmetry depends upon the initial value 
of phase, 9, of field x- According to the presented above arguments, 9 should 
be homogeneously distributed in the interval [0, 2tt]. Thus the asymmetry 
between x an d anti-% inside the bubbles may vary from zero to the maxi- 
mum value B x ~ m|x| 2 . The baryon asymmetry f3 depends upon the ratio 
of the energy density of x, Px ~ m2 |x| 2 > & t the moment of x decay to the 
background cosmological energy density. If x decayed when the background 
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relativistic matter was strongly red-shifted, the energy density inside the 
bubble would be dominated by the energy density of \ an d the ratio of the 
baryonic charge density to the entropy density inside the bubbles would be 
about 

*~<^ = gf- (18) 

This value may be much larger than unity, 
o.oi 

\%\ 

0.008 
0.006 
0.004 
0.002 



100000 200000 300000 400000 

t 

Figure 4: Evolution of |x| in time. The field rolls down toward the deeper min- 
imum, oscillates there following the evolution of the minimum, rolls back to the 
origin and starts to rotate around it. 




4 Distributions of bubbles over size and mass 

The process of the B-bubble formation predominantly took place when the 
effective mass of x- field, eq. (J7J), was smaller than the Hubble parameter (fTUl) . 
It roughly corresponds to 

M$ _ $l)2 <l^|£!. (19) 

3m 2 Pl 

During this period we consider x as a massless field and estimate how fast its 
quantum fluctuations rose. Here we follow ref. [IB], which in turn was based 
on the diffusion equation governing evolution of the quantum fluctuations, 
derived in paper [20] . 
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Qualitatively the picture looks as following. "Massless" field % "Brown- 
ian" moves from zero with the average displacement given by eq. (J2J). If X 
exceeds a certain value Xc, after which the motion of x away from the origin 
would be governed by its potential (JHJ), the motion turns into classical roll 
down to the second minimum X2 of this potential. The boundary value when 
the evolution of x turns from the quantum regime into the classical one can 
be found from the condition U"{Xi = H 2 , or approximately: 

< 2 pi 



Those x's which were able to reach this value would give birth to B-bubbles. 
Earlier formed B-bubbles would have considerably larger sizes because their 
creation took place at inflationary period and the ratio of the radii of two 
bubbles crated at t\ and ti would be: 

r 1 /r 2 =exp^ Hdtj . (21) 

Let us introduce the probability distribution V(x, t) of quantum fluctu- 
ations to have the value x — Xi + % at time t. We assume that it is 
normalized to unity: 



/+oo 
dxidx2V{x,t) = 1 . 
-oo 



:22) 



The probability of B-bubble formation per unit time and per unit volume 
can be parametrized as: 



dW 
dtdV 



J d X id X 2^V(x,t) , (23) 



where fi is a parameter with dimension of mass. It is natural to expect that 
H = H. There are simply no other parameter in the massless limit. At 
inflationary stage we may assume that H w const. The initial size of a B- 
bubble should be of the order of rt n ~ 1/ H. Hence the fraction of the total 
volume occupied by the bubbles is 

ft POO 

u(t)= / dt'H / dxidx 2 V(x,t') , (24) 

where we may formally take t in = — oo but in reality it is determined by the 
condition m e // < H. To obtain the above result we took into account that 
the bubbles formed at time t' exponentially expand together with the rest of 
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the universe and their radius rises as r(t) = H~ 1 a(t)/a(t'), where a(t) is the 
cosmological scale factor. So after formation their volume fraction remains 
constant. 

Correspondingly, we find that the volume distribution of B-bubbles, or 
better to say, their number density as a function of their size would have the 
form: 

POO 

N{r,t)~ d X id X 2V(x,t\r,t)) , (25) 
Jxc 

where t'(r,t) can be found from the equation: a(t)/a(t') = rH and H is the 
Hubble parameter during inflation. In particular, if we remain at inflationary 
stage, then a(t) = exp(Ht) and t'(r,t) = t — ln(rif ). This derivation is 
valid for constant H. In realistic situation H drops down with time. We will 
consider this case below. 

As is argued in ref. [TB] (see also this section below), for the quadratic 
potential U = m 2 e ^{t)\x\ 2 the probability distribution has the Gaussian form: 

V(x,t)=C(t)e W [-\x 2 \/2( X 2 )] , (26) 

where C(t) is determined by the normalization condition ([22)1 . The mean 
square width of the distribution is given by 

(X 2 ) = ^ jTcft'exp (-A jf dtV/,(*")) • (27) 

When $ is close to $1, the effective mass behaves as m\^{t) ~ m\ + mg(t — 
ti) 2 , where t\ is the moment when $ = $ x . This case was analyzed in ref. [16] . 
where it was shown that both the bubble distributions over length and mass 
have log-normal form. In particular, 

— = C M exp[- 1 \n 2 (M/M )] , (28) 

where Cm, 7, and M are some constant parameters, which can be expressed 
through the parameters of the model. For phenomenological applications we 
take Cm, 1, and M Q as almost free parameters with the values which are 
in the interval dictated by the "reasonable" values of the coupling constants 
and masses of x and $. Notice that a modification of this distribution by 
a power factor, M b , or, which is the same, by a log term in the exponent, 
exp(MnM) lead to the same log-normal form of the distribution (1281) with 
the corresponding change of the values of the parameters. 

The log-normal form of the bubble distribution was derived in ref. [E] 
under simlifying assumption of a constant Hubble parameter during inflation. 
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In many realistic situations H depends upon time and drops down practically 
to zero when inflation ends at t — t e : 



H(t)=H 1 h(r) = H 1 [l-h 1 (r-r 1 )]. 



(29) 



Here r = Hit and t\ is the moment when $ = $! and hi — l/(T e — Ti) <C 1. 
For such H(t) the size of the bubble created at t — t p rises as 



1(t,t p ) = / iri exp 



drh{r) 



(30) 



Correspondingly, the size of the bubble at the end of inflation will be 



K T e,T p ) = Z in exp 



(r e -r p ) [1 + hin - -/ii (r e + r p ) 



(31) 



From this equation we see that the bubbe size is related to the production 
moment as 



T p = hi 1 (l + h t Ti - ^2hilnlHi 



The distribution of the bubbles over their size can be written: 

2 



W(l) = Qexp 



7 [Te ~ n 



^2{\nl)/hi 



(32) 



(33) 



After some modifications this expression can be rewritten as 



W{1) = dexp 



- 7 (r e - ti) 1 1 



hi(r e - ti)' 



In/ 



(34) 



Now after simple algebra we find that the distribution over / again has the 
log-normal form: 



W{1) = Qexp 



2-i 



Mre-Tif 1-4/1 + 



\nl/l 
In In 



Qexp [-7 In 2 (l/k 



(35) 



where l = exp[hi(T e — Ti) 2 /2] and 7 = ^/h\{r e — Ti) 2 . 

The equation governing the evolution of the quantum fluctuations of a 
real scalar field was derived in ref. [20]. Following this derivation we can 
conclude that the probability distribution for a complex field \ in potential 
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U (%) in (quasi) De Sitter space-time with the Hubble parameter H satisfies 
practically the same equation 



~ 8^2 ^v 2 + ^77 2^ 



2.T) I A 

(36) 



dt 8tt 2 ^ 9x1 ' ^ a^fc ' '' 

fc=l,2 AA " fc=l,2 A 



P 



where Xi,2 are the real and imaginary parts of field, x — respectively. 
We decompose Xk into classical part x^ an d quantum fluctuations, 

Xk = xT + Vk ■ (37) 

At the initial stage, when $ was far from $ 1; the effective mass of x was large 
and quantum fluctuations were small near Xk = 0. In this case the potential 
governing their evolution can be approximated as the quadratic one: 

U q = l m2 ijViVj , (38) 



with the effective mass: 



2 _ d 2 U 

%3 dxidx. 



(0) 
Xi=X\ 



(39) 



Due to the last two terms in eq. (jHJ) the potential is not symmetric with 
respect to the phase rotation of x an d the matrix m^- may contain off-diagonal 
terms, m\i ^ 0. However by a suitable redefinition of x ~ > ex P(^0)x we can 
eliminate imaginary part of m and make diagonal. 

For a quadratic potential eq. (136!) can be solved in the form: 

V = A(t) exp(-a ij r]ir] j ) , (40) 

where symmetric 2x2 matrix and A is determined by the normal- 

ization condition (122j) : 

A(t) = -Vdeta . (41) 

7T 

After simple algebra, comparing terms of zeroth and second order in x%Xji 
we obtain in the matrix form: 

^r[ar l a] = -|^tr[a] + ^tr[m 2 ] , (42) 

lI a 2 + ^-am 2 . (43) 



2tt 2 3H 

There are three unknown functions, a>ij(t) and four equations (taken into 
account symmetry of matrix a). Still the system is not over-determined. If 
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we multiply the second equation by matrix a~ l and take trace, we obtain the 
first equation. 

As we mentioned above, the mass matrix m 2 can be made diagonal by 
phase rotation of x- Correspondingly matrix a becomes diagonal too and the 
equations for an and 0,22 components decouple. Introducing Kj = l/aj^ 1 (no 
summation), we obtain the linear equations: 

2 2 H 3 
^ + 3H m ^ = 2^- (44) 

These equations can be integrated by quadratures for arbitrary m 2 {t) and 
H(t) but it may be simpler to integrate them numerically directly (see Fig. 
ED- 



0.08 - 




5 10 15 20 25 



Figure 5: The variation of the mean square values function of time. All 

quantities are normalized to H\, i.e. t — ► Hit and kj —* kj/H 2 . This result is 
obtained with effective mass (j39|) using potential ([5]). 

This result allows to calculate the probability distribution for sufficiently 
small quantum part of \i when quadratic approximation is valid. As we 
have already mentioned, when \ reached the critical value Xc (ffi?]) its evo- 
lution would be governed by classical equation © with the potential which 
is dominated by log-quartic term A2 1 xl 4 l n \x\ 2 / a2 - These terms in the po- 
tential become essential rather early, still at the quantum regime. To get a 
qualitative understanding of the impact of this terms we proceed as follows. 
We approximate the potential by 

U^m 2 f \ X \ 2 -\2\x\\ (45) 
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where m 2 ,** = ml + Ai($ — $i) 2 . We neglect the last two aspherical terms 
in the potential (jHJ) because their role at this stage is not important. For 
a spherically symmetric potential U = U(\x\) the probability distribution 
depends only upon \x\ and we will look for the approximate solution in the 
form: 

V(x,t)=A(t)[l + b(t)\ X \ 4 } exp [-|xl7*(*)] , (46) 

assuming that terms containing b(t) with higher powers of x are small. 
The normalization condition now leads to: 

A = 7TK(1 + 2bK 2 ) ' ^ 



Substituting this expression into eq. ( 1361) and neglecting terms of higher order 
in b and sixth power of x we fhid for zero, second, and fourth powers of x 
respectively (only two out of these three equations are independent): 

k R3 + (ml ff +i\ 2 K) , (48) 



K 



2it 2 k v 1 3H 

2H 3 b . 8A 2 8m lff b 



It is easy to solve these equations numerically. As expected, k becomes a few 
percent larger at $ ~ $i, i.e. the additional negative quartic term stimulates 
the growth of fluctuations. 



5 Density perturbations created by B-bubbles 
and their evolution 

As it was mentioned in ref. [IB] . B-bubbles could create large isocurvature 
density perturbations at relatively small scales. The spectrum of angular 
fluctuations of cosmic microwave background radiation (CMBR) excludes 
noticeable isocurvature perturbations with wave lengths larger than approx- 
imately 10 Mpc [2SH2Z], but at smaller scales they are allowed and may even 
be very large. 

In the considered model there are two mechanisms of the density per- 
turbations, which were operative at different cosmological epochs. Inside 
B-bubbles the baryon asymmetry, /3, may be much larger than the observed 
value f3 ~ 6 • 10~ 10 . Such inhomogeneous distribution of the baryon number 
would not lead to noticeable density perturbations till quarks remained ultra- 
relativistic, i.e. till the QCD phase transition. In this way relatively large 
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objects with high density contrast could be created. This was considered in 
ref. [16] , where stellar and larger mass objects were studied. There is another 
possibility to perturb the cosmological energy density disregarded in ref. [16] . 
This mechanism operated at an earlier cosmological period and might lead to 
perturbations with much smaller wave lengths and masses, M <C M , than 
those created at the QCD phase transition. Here we consider in detail both 
possibilities. 

Let us first discuss the early stage when the B-bubbles were created. 
Initially, at inflationary stage, the density contrast between low and high 
X-regions was negligible, first, because the energy density of % is assumed 
to be strongly subdominant with respect to the inflaton one. Second, the 
bubble formation did not lead to a density contrast at the moment of for- 
mation simply because of energy conservation. The energy density contrast 
Sp between the bubbles and the rest of the universe remained small till the 
end of inflation. After inflation is over and the inflaton decayed into (pre- 
dominantly) relativistic particles, the situation might change, if the life-time 
of x is larger than the time of the inflaton decay. Outside the bubbles there 
would be only relativistic matter, while inside the bubbles there would be a 
mixture of nonrelativistic and relativistic matter. Indeed, in the course of 
the subsequent cosmological evolution the equation of state of x would be 
different for small x near zero and for large x staying in the second minimum 
of U (x, $) ■ At this stage the energy density of x would behave as vacuum- 
like energy. When this minimum disappeared, x started to oscillate around 
zero in essentially quadratic potential. So its equation of state became that 
of nonrelativistic matter and its energy density was red-shifted as the scale 
factor cube. 

To illustrate the evolution of the density perturbations in the considered 
model let us assume that there are two large (larger than horizon) pieces of 
the universe, the bulk with the relativistic equation of state, p r = p r /3, and 
a bubble with a large value of x with a mixture of relativistic and a small 
fraction of non- relativistic matter, p™ r = ep™. The energy density inside the 
bubble evolves as: 



„in „in „m 



p insid e = ^ I + ^ = ^ T (l + ex) , (50) 



X 4 X 3 X 4 



where x = a/a>i n is the ratio of the running cosmological scale factor to its 
initial value. 

We can find the dependence of the scale factor on time solving the Fried- 
mann equation: 



x 



2 



3m? pi V x A x 3 



' N, ' t + ^ nr ) (51) 
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Technically simpler is to solve this equation in terms of conformal time, 
dr = dt/a(t): 

a = a eq (2t] + 7] 2 ) , (52) 

where a eq = a in /e is the value of the scale factor when the energy densities 
of relativistic and nonrelativistic matter inside the bubble become equal and 
rj = Tai n /(Ati n e). At matter radiation equilibrium rj ~ 1. 

In what follows we will use the results presented in book [28] where a very 
similar case of evolution of density perturbations in the mixed relativistic and 
nonrelativistic matter is considered. The gravitational potential $ satisfies 
eq. (7.49) of this book with the right hand side which is non-zero because 
of the difference of equations of state inside and outside the bubble. The 
solution for $ with the initial condition $(77 = 0) = has the form: 

, pf + eri+W SS 
* = 5 (, + 2)' ~S ' (53) 



where 

5S 3 5p r 5p„ 
S 4 p r p m 



(54) 



Notice that the entropy perturbations are not small becasue Sp nr /p nr = 1. 

The density perturbations are related to the potential through the equa- 
tion (7.47) of book 



A$ - m ($' + H<$>) = AnG N a 2 5p, (55) 

where the Laplacian term in the l.h.s. can be neglected for large wave length 
of perturbations. Thus Sp/p ~ 2$. 

From this equation we find that asymptotically at MD stage, i.e. for 
r] ^> 1, the density perturbations tend to a constant value 



Sp 
P 



2 SS 2 fK^ 

= -7^~7- (56) 

MD ° ° ° 



For small 7], i.e. at RD stage the perturbations rise linerly with conformal 
time: 



Sp 
P 



= -I f - f • (57) 

RD Z D Z 



At some moment field \ should decay into relativistic matter and after 
that both inside and outside the bubble equations of state would be identical, 
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those of relativistic matter, p = p/3. After decay the developed density 
perturbations would stay constant, see eq. (7.59) of ref. [28J. In the instant 
decay approximation the initial value of the constant potential and density 
contrast can be read off equations ( 1561) and ( 1571) . To be more accurate we 
can proceed as follows. We have to solve the system of the following two 
equations, assuming that all nonrelativistic matter is created by x-field: 

p x + 3H Px + T xPx = , (58) 
p r + 4Hp r - T xPx = , (59) 

where T x = l/r x is decay rate of field x m t° light particles. 
The solution of equation ( 1581) : 

in 

p x = ^exp(-r x t). (60) 

Equation (|59[) would be simplified in terms of the new function p r = p* n £ / x 4 
with the initial value £ = 1 at t — 0, or x — 1. The resulting equation takes 
the form: 

C = er x xe- r *', (61) 

where x = x(t) is determined by eq. (|5"Ti) but instead of p given by eq. (1ST)]) 
we have to substitute there the solutions of eqs. (I58|59p . A solution to this 
system of equations can be done numerically. 

We see that the density contrast for superhorizon perturbations initially 
grows up but asymptotically tends to a constant value. When the bubble 
size becomes smaller than cosmological horizon, the further evolution of the 
bubble depends on the magnitude of the density contrast at the moment of 
the horizon crossing, A^. If Ah ~ 1, such bubbles would form primordial 
black holes (PBH). For a review of PBH creation see e.g. refs. [221 ED]- 

We should remember however, that the presented above estimates are true 
only for small perturbations when first order expansion in perturbed quanti- 
ties is accurate. As we saw, the estimated perturbations are of the order of 
unity and so we can trust the calculations only by an order of magnitude. 

The mass inside the horizon is 

M h ~ m 2 Pl t . (62) 

So the discussed mechanism allows for formation of PBH with the masses of 
the order of M BH ~ 4 • 10 38 (r x /sec)g, if the size of the bubble was inflated 
up to horizon at that moment. An oversimplified estimate of the life-time of 
X-boson as r x ~ (g 2 m x /ATr)~ 1 , where g is the Yukawa coupling constant of \ 
to fermions, gives for g 2 = 0.01 and m x = 10 8 GeV: Mbh ~ 10 9 g. Such light 
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PBHs should evaporate in about one second and would not be important in 
the present day universe. On the other hand, the life-time estimated above is 
certainly too short. In supersymmetric theory with conserved 7?.-parity the 
lightest supersymmetric particle (LSP) must be stable. So x should either 
be stable or decay into a massive LSP. For simplicity we will not dwell on the 
problems related to color conservation and formation of x condensate. This 
discussion can be found e.g. in the reviews [311 E2]- We will assume that 
7?.-parity is broken and x predominantly decays into ordinary particles but 
with baryonic charge conservation. If this is realized, the life-time may be 
much smaller than the estimate presented above and PBHs with much larger 
masses could be created. Since 7£-parity violating effects are only bounded 
from above, the life-time is allowed to be arbitrary large. Another effect 
which might result in larger PBH masses is a delayed evaporation of the x~ 
condensate found in refs. (33J, [33]. On the other hand, it is not obligatory 
to live in a SUSY world but instead we can postulate ad hoc an existence 
of a new %-field with the necessary properties and in this case we would 
be free from the SUSY restrictions. So, one way or another, the considered 
mechanism can lead to formation of PBH which live long enough to be present 
in the contemporary universe. 

If the density contrast is not sufficient for creation of PBH, the bubble 
would end up as a more dense piece of relativsitic matter inside the relativistic 
cosmological backgound. The density contrast 5 inside horizon is known to 
remains constant or to be more precise, may rise logarithmically (see e.g. 
[35l l28l 36J), while outside horizon it drops down as 1/t, as is mentioned 
above. The bubbles might be dissolved by diffusion but in relativstic regime 
the mean free path is about lf ree ~ l/(a 2 T) and thus remains small in 
comparizon with the bubble size which rises with the cosmological scale factor 
a ~ 1/T. So if the bubble size was initially larger than lf ree it remains such 
during all relativistic stage. 

The second period when the perturbation in baryonic number density 
could be transformed into the energy density perturbation, happened after 
the QCD phase transition (p.t.) from deconfinement to confinement phase. 
After this p.t. practically massless quarks turned into heavy baryons and the 
baryonic number fluctuations would be transformed into the energy density 
perturbations. If the baryon asymmerty inside the bubble has the value (3 b 
the density contrast after the QCD p.t. would be 

A = — = -. - 7 X P ^ A w 0.07fi B -T^r , 63 
p (tt 2 /30)^T 4 hu T ' v ; 

where temperature T should be smaller than the temperature of the QCD 
p.t., Tqcd- The latter is rather poorly known and according to different 
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estimates it varies between 100-200 MeV. The quoted results, however, were 
obtained for zero or negligible baryonic chemical potential, \xb- In our case 
Hb could be large, but not large enough to shift the value of Tqcd for more 
then 2% [33 EH]. 

The bubble with A (1631) larger than unity at horizon crossing would form 
PBH but now (in contrast to the case consided above) with very large masses, 
starting from a few solar masses to possibly millions of M , if the bubble 
size is so large that it enters horizon at, say, t — 10 — 100 s, see eq. ( J62l) with 
t > 10~ 5 s. 

If at horizon crossing A# < 1, the bubble evoluton is detemermined by 
the relation between its size, Ib and the Jeans wave length: 



irm 



pi 



P 



1/2 



(64) 



where p is the background energy density and c s is the speed of sound. For 
relativistic plasma c s = 1/ a/3 and the Jeans wave length is close to horizon 

lh = 2t. 

The Jeans mass is defined as 

Mj = —p\j 3 . (65) 

As is well known, the regions with M > Mj or with A > Aj would decouple 
from the cosmological expansion and form gravitationally binded systems. 
In the regions with smaller mass gravitational instability is not developed 
and the perturbations oscillate. 

In the usual cosmological plasma the Jeans mass and wave length are very 
large before hydrogen recombination because of very small baryon asymmetry 
or a very large entropy per baryon, see e.g. book [39]. After e + e~ annihila- 
tion, i.e. at T < 0.1 MeV the baryonic Jeans mass was about 1O 4 M and rose 
up to 1O 18 M till hydrogen recombination. After recombination it sharply 
drops down to 1O 6 M and continue to decrease as a(t)~ 3 / 2 . 

In our case the entropy per baryon is not huge and the picture is quite 
different. The energy density inside B-bubble, assuming that baryons are 
nonrelativstic, is 

p B = n B m p + (7r 2 /30)g*T 4 , (66) 

where % is the baryon nuber density and g* is the efective number of the 
relativistic species. For photons g* = 2 and for equilibrium photons and e ± 
pairs g* = 5.5. 
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The pressure density is equal to: 



P B = (tt 2 /90)^T 4 + n B T . (67) 

For adiabatic perturbations the ratio of the baryon number density to 
the photon number density remains constant if entropy release by e + e~- 
annihilation is neglected. Disregarding this and similar efects we obtain by 
an order of magnitude: 

C ' = ItT^- < 68 > 

Correspondingly the Jeans mass at the QCD phase transition would be about 
10M Q if /3 = 1. Assuming that the B-bubble is dominated by nonrelativsitc 
matter we find that with further cosmological expansion Mj would drop as 
a 3//2 because the energy density of baryons decreases as a~ 3 and the tempera- 
ture of mixture of nonrelativistic baryons and photons decreases as T ~ 1 /a 2 . 
For smaller f3 the Jeans mass can be significantly larger. 

Thus we see that after the QCD p.t. PBH with masses starting from a 
few solar masses and much higher can be created as well as gravitationally 
bound objects with masses from a fraction of the solar mass up to thousands 
or more solar masses could emerge. These early formed objects would be 
made either from matter or antimatter with equal probability. 



6 Observational bounds and implications 

There are two features of the considered here model which may make it 
cosmologically and astrophysically interesting. 

First, according to the presented scenario all cosmological dark matter 
can be made out of baryons and antibaryons in the form of compact stellar 
type objects or primordial black holefl These PBHs can be rather light, 
e.g 10 15 — 10 20 g or very heavy, up to millions solar masses. The number 
density of light PBHs is restricted by the condition that their evaporation 
should not create too strong electromagnetic radiation to exceed the observed 
one. The analysis, summarized in ref. [291 |4"T] excludes PBHs in a large 
interval of masses as the dominant part of cosmological DM. However, the 
conclusion of ref. [JT] is valid only for a narrow spectrum of PBH masses, 
when all PBHs have practically equal masses. In our case the mass spectrum 
can be sufficiently wide to allow for light PBHs to be the dominant dark 

1 A different model of cosmological dark matter made of baryons in the form of astro- 
nomically small quark lumps, ~ 10 7 g is suggested in ref. [30] 
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matter (DM) particles and to avoid condradiction with the observed cosmic 
electromagnetic radiation. 

A possible manifestation of PBH with masses around 10 16 — 10 17 g may be 
0.511 MeV line observed in the Galactic bulge pEl H3J Hfl SHI EE S3 HE] and 
probably from the halo [10]. As argued in ref. [50] such PBHs may explain 
the observed flux and be abundant enough to make all cosmological DM. As 
was observed recently the source of this line is asymmetric with respect to 
the galactic center and is correlated with the observed X-ray binaries [5*T] . 
It make probable that the origin of the annihilated positrons are these low 
mass X-ray binaries, though the estimated intensity of the line is about a 
half of the observed. On the other hand, the distribution of PBHs may be 
also asymmetric and correlated with the distribution of the binaries. 

If the central mass in the distribution (|28p is close to the solar mass, we 
would expect quite heavy cosmological dark matter "particles" made out of 
heavy PBH and dead or low luminocity stars with masses of a fraction of 
solar mass and may be much larger. An analysis of the evolution of such 
stellar like objects formed in the early universe will be made elsewhere. As is 
argued in ref. [TB], an interesting feature of this scenario is that on the high 
mass tail of the distribution (128j) there may be superheavy PBH with masses 
of millions solar masses. We can choose the parameters of the distibuiton 
such that the number of the superheavy BH would be equal to the number of 
the observed large galaxies. So the model explains an early quasar creation 
which may seed galaxy formation, as it was recently argued [52J. To make a 
better qualitative statement one should study the amount of matter accretion 
during life time of these superheavy BHs from the moment of their creation 
when the universe was about 1 sec old till practically present time. 

The observed in microlensing searches stellar mass low luminocity objects 
(MACHOS) could be such solar mass PBH or dead very early stars and 
antistars. According to ref [53], such objects may make at most 40% of the 
total galactic DM. However, again this bound is true only for narrow mass 
spectrum of MACHOS. If the latter is sufficiently wide PBH can be safely 
put below the upper bound curve and make 100% of galactic DM. (A little 
here, a little there makes unity.) 

An argument against very heavy DM particles is presented in ref. [54J, 
according to which the objects with masses above 1O~ 2 M may constitute 
not more than 88% of cosmological DM at 95% confidence level. So we should 
either allow for 12% to be in the form of relatively light PBHs or planet-like 
objects with masses below 2- 10 31 g or hope that 95% CL is not strong enough 
to exclude very heavy PBHs. 

The second possibly testable feature of the model is a prediciton of no- 
ticeable amount af antimatter, in the form of compact stellar-like objects. 
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Out of VLdm — 0.25 a half may consist of antimatter. So we have almost 
charge symmetric universe with the amount of matter fig ~ 0.16 and that 
of antimatter Q§ « 0.12. Of course the bulk of baryonic and antibaryonic 
matter can be inside PBHs and thus unboservable (except for its gravity), 
however, it is not excluded that a noticeable amlount of antimatter is in 
the form of evolved anti-stars or antimatter clouds. There are quite strong 
bounds on the amount of cosmic antimatter in the usual form. For example, 
the nearest antigalaxy mast be at least at the distance of 10 Mpc [55J. In 
baryosymmeteric cosmology the nearest domain of antimatter should be at 
the distance of (a few) Gygaparsecs [13] . Bounds on the amount of galactic 
antimatter in "normal" form in baryo-asymmetric cosmology are discussed 
in ref. JHUIEZIEH]- However, in our case these bounds are not applicable, see 
ref. [2], and one may expect that roughly a half of the galactic mass is made 
of antimatter in the form of compact low luminocity stars or PBH situated 
in the disk and in the halo. Possible antistars in the galactic Bulge may be 
sources of the observed positrons [2]. 

In the regions with high (anti)baryonic number density primordial nu- 
cleosynthesis (BBN) would lead to quite diferent amount of light elements. 
There would be practically no deuterium, somewhat larger amount of A He 
and 7 Li, and anomallously larger abundances of heavier elements which prac- 
tically are not produced in the standard BBN j59[ EH [61]. This feature can 
indicate that a cloud or a stellar type objects with anomalous light element 
abundances may consist of antimatter. 

According to ref. [62l [63] the moment of recombination is sensitive to the 
spatial distribution of baryonic matter and it might be helpful for restriction 
of the model parameters or its confirmation. 

For very simple form of the interaction between x an d the inflaton field 
Q the log-normal mass spectrum of BBs has maximum at some M which 
can be either close to normal stellar mass or much smaller, say, 10 15 — 10 20 g. 
However, with a minor modification of the discussed model we can obtain a 
multi- maxima spectrum with several values of Mq. To achieve that we need 
to modify the interaction of x with the inflaton as 

U int = \ k U){^-^ j f/m^- 1 \ (69) 

Such an interaction may emerge from non-renormalizble Planck scale physics. 

7 Conclusions 

We reconsidered and further developed the scenario of ref. [16] and argued 
that a simple coupling of the Affleck-Dine scalar baryon to the inflaton may 
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lead to an efficient creation of PBHs either with relatively small masses, 
Mbh = 10 14 — 10 20 g, or much larger, M BH > M Q . The predicted log- 
normal mass spectrum of PBH ( |28l) very much differ from the power law 
one which could be created by other mechanisms considered in the existing 
literature. This difference is crucial for relaxing the observational bounds 
on PBH abundance. The suggested mechanism may explain formation of all 
observed heavy black holes from millions to thousands or tens solar masses. 
There is no unique and compelling explanation in the existing literature. 
There is an interesting possibility that all cosmological DM is made of PBH 
or compact stellar mass objects (early formed and now dead stars) with log- 
normal mass spectrum. In particular, the observed MACHOS may be solar 
mass PBHs, or dead (anti) stars. 

In the simplest version the mechanism considered here leads to signif- 
icant creation of astronomical objects made of antimatter. Such objects 
should have anomalay abundances of (primordially formed) light elements. 
The search for cosmic antimatter |64J, if successed would be crucial for the 
conformation of the scenarion. 
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